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1. INTRODUCTION 

The environmental and climate concerns caused by generators using fossil fuels are some of the 
reasons for an urgent need of the world to consider introducing alternative energies from solar, wind, 
geothermal, and others into its existing power production [1]. Solar photovoltaic is one of the widely used 
renewable energy technologies in the world. But nonlinear loads and converters commonly available in PV 
systems cause harmonic currents. Harmonics is a critically important concept in electronics. Harmonics 
affect the quality of AC electricity delivered to homes and facilities and the performance of equipment that 
uses this electricity. Harmonics can increase energy costs and reduce the lifespan of the hardware. In some 
cases, harmonics can overheat electrical conductors, creating a fire risk. Additionally, the control of the 
voltage source inverter (VSI) is highly affected by the presence of harmonic currents which contribute to the 
poor power quality in the system. The current distortion level is commonly indicated by the total harmonic 
distortion (THD) given by [2]: 


o Th 
THD; = (1) 


n=2 lL 
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Where: 
lL = fundamental current. 
[,= current at n™ harmonic 


Several techniques have been applied to damp harmonic currents in a grid-connected PV system. 
The methods use filters to connect the VSI with the grid and controllers to control the VSI. The filters applied 
are such as L, LC, LCL filters which interface the VSI with the grid. The conventional proportional-integral 
controller, fuzzy logic control, optimal control, and robust control have been widely used to control the VSI. 
Most of the optimal and robust control methods have been designed based on the dynamic model of the LCL 
filter. Additionally, low pass and notch filters have been part of these control strategies. However, there is a 
lack of h-infinity methods for harmonic current distortion. Thus, this research aims at bridging this gap in the 
existing literature. 

From the literature review of the existing techniques, some presented below, this paper presents a 
new current compensation control method that focuses on the reduction of the grid current THD. The 
contribution of this research work which has not been done before is an improvement in the THD of the PV 
system by applying robust loop shaping control for the voltage source inverter (VSI), filtering actions of an 
LCL filter, and active damping using notch filters. The results obtained after simulations show the robustness 
of the proposed method in damping out harmonic currents. 

L. El Iysaouy, et al., in [3] compare the performance of CL and LCL filters in THD with both filters 
manage to reduce the THD to less than 5% fulfilling the requirements of the IEEE 1574-2014. Authors in [4] 
added more functions to the proportional resonance (PR) current controller to minimize THD of the 3rd 
harmonic from 0.45% to 0.1%, the 5th harmonic from 0.6% to 0.25%, and the 7th harmonic from 0.43% to 
0.4%. In [5] a control strategy based on parallel-connected LCL-type inverters is used for harmonic voltage 
distortion. Passive damping is added to the LCL filter in [6] for improving the power quality in a solar 
photovoltaic system. Not only the grid voltage THD is reduced from 0.59% to 0.08% but also the grid current 
THD has dropped from 10.71% to 1.17% and the inverter current THD is reduced from 10.70 to 4.9%. 
Electromagnetic interference (EMI) filter and LCL filter are combined to form a hybrid technique named 
EMI-LCL filter for harmonic and EMI noise suppression for single-phase SiC-MOSFET grid-connected 
inverter [7]. 

W. A. A. Saleh, et al., in [8] apply the nearest level control (NLC) method to a 13-level transistor- 
clamped H-bridge (TCHB) inverter to reduce the harmonic content because this control structure consists of 
two TCHB cells supplied with two asymmetrical DC input sources reduces the number of electronic 
components in the inverter. The optimal setting of the inverter’s filter elements against irradiance is 
performed in [2] to eliminate harmonics at low irradiance levels. The test is conducted under faulty 
conditions. 

M. A. Razak, et al., in [9] manage to minimize 5th, 7th, 11th, and 13th harmonics caused by 
semiconductor devices and converters in a PV system by using a passive filter. A Hammerstein adaptive 
filter is proposed in [10]. The method is capable of extracting the positive sequence components of the grid 
voltage during imbalance. This parameter is used to estimate the unit templates and the proposed control 
strategy is left for load current fundamental components extraction. The control strategy limits the THD to 
the IEEE-519 standards. 


2. PV SYSTEM WITH AN L FILTER 
The grid-connected PV system which includes the PV source, the DC/DC converter, the VSI, and an 
L filter is shown in Figure 1 [10]. 











PV Module DC/DC Converter DC/AC Converter Filter Grid 


Figure 1. A grid-connected PV system with an L filter 
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2.1. PI control 
The reference gate signal voltages in dq transformation to control the VSI in Figure 1 by using the 


PI control strategy are given by (2) and (3) and explained in detail in [11]-[13]. 
Viaref (t) = Vga — LWsigg + Kpea + K; J eqdt (2) 
Viqref (t) = Vgq + LWsiag + Kpeg + K; J eqdt (3) 
Where ea = iaref — iga and eg = igref — İgq 


2.2. Fuzzy logic control 
The fuzzification, fuzzy rule base, and defuzzification steps are shown in Figure 2. 


ligref + 












Fuzzification Difuzzification 





Figure 2. The fuzzy logic control structure 


Where: 

e(t)= Error 

é(t) = Derivative of the error 
u(t)= Output signal 


The design of the "if -and —then"rules is based on qualitative knowledge, deduced from extensive 


simulation tests using the trial-and-error method. 

In this scheme, the PlI-controllers used in the previous method in the reference voltage calculation 
process are replaced by fuzzy logic controllers as shown in Figure 3 [14], [15]. The errors e(t) (obtained 
from the subtraction of the measured grid currents from their reference currents) and their derivatives é(t)are 
fed to their corresponding fuzzy controllers which output voltage signals summed up with the measured grid 
voltages and cross-coupling components to obtain dq -axis reference voltages. Then, these dq -axis voltages 
are converted into the abc reference frame. 





Figure 3. FLC of current loops 
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The fuzzy rules are as shown in Table 1 and defined as follows: negative big for NB, negative 
medium for NM, negative small for NS, negative very small for NVS, zero for ZE, positive very small for 
PVS, positive small for PS, positive medium for PM, and positive big for PB. 


Table 1. Fuzzy logic rules 
e(t) 
e(t) NB NM NS ZE NS NM NB 
NB NM NM NM NS NS NS ZE 
NM NM NM NNS ZE NS ZE PM 
NS NB NM NS NS ZE PB PS 
ZE NM NN NM ZE PB PS PM 
PS NB NB ZE PS PS PM PB 
PM NM ZE PM PM PS PM PB 
PB ZE PB PB PM PM PB PB 


3. PV SYSTEM WITH AN LCL FILTER 

In this configuration, instead of an L filter used in Figure 1, an LCL filter interfaces the PV system 
with the grid as shown in Figure 4 [16], [17]. 

The filter in Figure 4 is represented by its equivalent circuit in Figure 5 [18] and its transfer function 
block diagram in Figure 6, V; which denotes the inverter voltage, V, denotes the grid voltage, w, is the grid 
voltage angular frequency, Cand v, is the filter capacitance and capacitor voltage respectively, i; is the 
inverter current, i, is the grid current, R;, Rg, Li. Lg is the filter resistances and inductances, respectively. The 
current is flowing from the inverter to the grid. 


oA alga vca Bla Ala 
ol B bla Vçb ab | B lo 
ac clo Vcc olc = Clo 
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Figure 4. An LCL filter is connected to the grid 


The Bode plot of the filter in Figure 7 shows that the filter generates a large overshoot at the 
resonance frequency. This overshoot results in the amplification of the harmonic amplitude at the resonance 
frequency, thus increasing the content of higher-order harmonics in the grid current. Therefore, damping 
resistance connected in series with the capacitor can be used to decrease the overshoot and resonance in the 
filter while maintaining the merits of the filter [16], [19]. This method is called passive damping. 

But, the presence of the damping resistance contributes to the increase in losses because the current 
Iç flows through it as explained in [20]. The losses increase as Rg increase. Additionally, the performance at 
steady-state and during transients is poor. To overcome this drawback active damping based on notch-filter is 
introduced and is used in this work instead of the damping resistance. 





Figure 5. LCL filter equivalent circuit Figure 6. LCL filter plant model 
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Signals in a narrow band around a given frequency are attenuated by a notch filter [17]. It can be 
tuned to reject frequencies near an LCL filter's resonance frequency, removing the resonant peak. Generally, 
the notch filter's transfer function is given by, 


_ s*4+2&wnst+w4 
Gy(s) 524.2% Wnstw2 (4) 


Where: 
W,= Natural frequency 
€,,é,= Damping factors 


The Bode plot of the LCL filter with and without damping is shown in Figure 8 [21] with the notch 
filter transfer function being. 


s7+25.945+5.286X10°8 (5) 
s*+2530s+5.286X10'8 


Gy(s) = 
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Figure 7. Bode plot of an LCL filter Figure 8. Bode plot of the undamped and damped 
LCL filter 


3.1. H-infinity control 
The block diagram representation of h-infinity optimization is shown in Figure 9 [22]. 








Figure 9. The standard Hoo configuration 


Where: 
w = Disturbance input 
u = Control input 
Z = Error output to be minimized 
y = Output of the plant 


The method seeks for the controller K to minimize the error zZ. 
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IITwz(P, K) Il (6) 


Subject to K stabilizes P. 

The objective function is the Hoo-norm of the closed-loop performance T,,_,,(P, K). This reduces the 
Hoo-norm of the transfer function from w to z less than a pre-determined positive number y to achieve system 
stability. The Hoo-norm can be mathematically expressed as; 


ITw2(s)llo < y (7) 


3.1.1. H-infinity controller implementation 

An LCL filter plant model shown in Figure 6 is used to synthesis the controller K. The synthesis 1s 
carried out by the command “hinfsyn” available in the robust control toolbox of MATLAB and implemented 
in SIMULINK as shown in Figure 10 [23]. 





Figure 10. LCL filter plant model with the Hoo current controller and notch filter 


3.2. Robust h-infinity loop shaping control 
The configuration of the loop shaping control with the plant P, controller K, input r, output y, and 
disturbance d and noise n is shown in Figure 11 [24]. 


l a pfe P -e a 4 
_ | n 








Figure 11. Loop shaping control block diagram 


The open-loop transfer function of the system is given by: 

L=PK (8) 
The output Yy can be written as, 

y = Pid +PK(r— y-n) (9) 
The error is, 

e=r-y-n (10) 
Multiplying both sides of (9) by identity matrix and solve for the output y gives, 

y = (I + PK) *PKr + (I + PK)*P,d — (I + PK)"*PKn (11) 
The sensitivity of the system can be given by, 


S= (I+ PK)! (12) 
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Complementary sensitivity can be given by, 
T = (I + PK)*PK (13) 
Substituting (8) into (12) and (13), Sensitivity and Complementary sensitivity can be written as, 
S=(+L)" (14) 
T=(1+L)1L (15) 
Let the addition of sensitivity and complementary sensitivity equal to the identity, 
S+T=I1 (16) 
Therefore, the error can be written as, 
e=Sr—SP,d+Tn (17) 
The frequency response of the loop shaping controller is shown in Figure 12 [24] and must satisfy 


the following conditions: 1) strong amplification in the low-frequency region, and 11) weak amplification in 
the high-frequency region. 


(log sa 
arale) aa L (jo) upper bound for 
ʻ high frequencies 
na : ie ii 
0 dB NN m 
lower bound for a 


low frequencies 





Figure 12. Robust bound of the loop shape 


The above frequency response ensures: 1) stability improvementii, 11) uncertainty compensation, 111) 
disturbance rejection, and iv) noise attenuation. 
The requirement for y tracking r is 


L 


yar El >>1 (18) 


14+L 


3.2.1. Loop shaping controller synthesis 
An LCL filter plant model shown in Figure 6 is used to synthesis the controller. The synthesis is 
carried out by the command "loopsyn" available in the robust control toolbox of MATLAB and implemented 


in SIMULINK as shown in Figure 13 [23]. 





Figure 13. LCL filter plant model with the loop shaping current controller, and notch filter. 
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The steps in the design of the controller are: 

1) Specify target loop-shape transfer function whose performance and robustness bounds are as in 
Figure 12. 

2) Use the target loop-shape to shape the nominal plant (Figure 6) to make its singular values match 
that of the target loop shape. 

3) Compute the loop-shaping controller ‘kK’ and the number 'y' that makes the singular value plot of the 
shaped loop matches the target loop shape where ‘y’ is in the range of 1 < y < 3 and indicates the 
accuracy to which the optimal loop shape matches the desired loop shape. 


The chosen desired loop shape was Gg(s) = ` and the obtained number y = 1.4163 which is within 


the specified range. The singular value plot of the plant with the controller is shown in Figure 14. The plot 
shows that the shaped plant optimally tracks the desired loop shape. Moreover, it can be seen from the upper 
half of the plot that the singular value plot of the open-loop gain best fits with the reciprocal of the singular 
value plot of the sensitivity (S) function, and in the lower half (below O dB line) the singular value plot of the 
complementary sensitivity (T) function matches that of the open-loop gain. All these ensure good 
performance, reference tracking, and disturbance rejection. 
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Figure 14. Singular values plot of the plant with the controller. 


4. RESULTS AND DISCUSSION 

The validity of the system is checked with the controllers implemented in MATLAB/SIMULINK 
for computer simulation to simulate the final system for different circumstances to analyze the system 
behavior. The use of a sim power system, Mamdani fuzzy inference system, and robust control toolbox 
provided the flexibility in designing the high rating system for actual implementation. The simulations were 
run on a PC having Intel Core 15, a 3.20 GHz processor with 8 GB RAM. The sample time is taken 20us. 
The effectiveness of the controllers is checked by considering two cases. In each case, an analysis is carried 
out and results are recorded for stability analysis. The PV system parameters are shown in Table 2. 

Case 1: The system is simulated for 2s under full load conditions which includes linear and 
nonlinear loads. 

Case 2: This case examines how the switching effects contribute to the increase in harmonics in the 
system and the response of the control schemes against switching transients. The system is simulated initially 
with linear load and then nonlinear load switched on and removed later. 


Table 2. Parameters of the PV system 


Parameter Value 
Grid voltage V,_, 380 V 
DC bus voltage vac 650 V 
Grid frequency f 50 Hz 
Inverter-side resistance R; 0.01 2 
Inverter-side inductance L; 0.7 mH 
Grid-side resistance Rg 0.001 2 
Grid-side inductance Lg 0.04 mH 
Filter capacitance C 50 uF 
Damping resistance Ra 0.0001 2 Q 
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4.1. Results-constant linear and nonlinear load 

The system is simulated for 2s for constant linear and nonlinear load. Simulations were recorded 
between T=1.9s to T=2s for all control schemes. During PI control the THD of the grid current was 7.85% as 
shown by the FFT analysis plot in Figure 15 and its corresponding grid current waveform in Figure 16. After 
the application of the fuzzy logic control, the distortion of the grid current dropped to 2.72% as shown in 
Figure 17 and its grid current in Figure 18 showing significant improvement. H infinity control managed to 
further reduce the THD to 1.52% as shown in Figure 19. At this point, the grid current ripples remain a few 
(see Figure 20). The application of the robust h-infinity loop shaping control removed all the distortions of 
the grid current leaving the waveform almost smooth as can be seen in Figure 22 with a reduction in THD 
dropping to 0.79 (Figure 21). The simulation waveforms show that the system with robust h infinity loop 
shaping controller’s performance is better in terms of dynamics and grid current smoothing. Moreover, the 
loop shaping-based controller is faster and reduces ripples in grid current and hence lesser deviation from the 
reference value. 
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Figure 15. Harmonic spectrum-PI control 
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4.2. Results-constant linear load and switching of nonlinear load 

The system is simulated from 0.95s to 1.1s for constant linear load and switching of nonlinear load 
between 1.00s and 1.05s. The simulation results for all controllers are shown in Figures 23 to 30. The 
nonlinear load is switched on at 1.00s and switched off at 1.05s for obtaining the switching transient response 
of the system. The simulation waveforms show that the system with a robust loop shaping controller’s 
performance is the best in terms of dynamics and grid current smoothing. Figures 23, 25, 27, and 29, show 
the performance of the simulated harmonics of the system during the switching of the nonlinear load. These 
harmonic spectrums show that the harmonics increase due to the switching transients of the nonlinear load. In 
this case, the transients have increased the harmonic contents of the grid current approximately 1.45 times the 
values recorded in the first case. With PI control the THD is 11.67% from 7.85% recorded in the first test 
(see Figure 23). This is verified by the system response when using the PI controller. The scheme is the 
poorest as the switching of the nonlinear load at 1.00s caused the shooting of the current to about 5.2A for the 
blue phase and about -5A for the red phase (see Figure 24). 

This was followed by approximately near steady-state values but with big ripples in the current 
waveform for all three phases. When the load was removed from the system at 1.05s, no switching transients 
occurred. The waveforms became smooth but slightly unstable. Fuzzy logic control managed to reach the 
allowable distortion margin set by IEEE-519 standards (THD of not more than 5%) because with this scheme 
the THD is 3.98% from 2.72% of the previous test (see Figure 25). The grid current waveform’s red phase is 
-4.5A and blue phase shot to 4.5A from its reference value due to the switching of the nonlinear load and 
slight variations of the green phases can be seen in Figure 26 followed by dropping of the current to the 
steady-state value with some ripples lower than those of the PI control. When the load was switched off at 
1.05s the grid current gained stability and the waveforms became smooth again for all three phases. H- 
infinity controller and robust loop shaping give 2.22% and 1.14% of THD respectively despite the switching 
effect of the nonlinear load (see Figures 27 and 29) although higher than the values obtained in the first case 
which were 1.52% for h-infinity control and 0.79% for the loop shaping control. The grid currents shown in 
Figures 28 and 30 prove the results of the harmonic spectrums in Figures 27 and 29. The grid current when 
the system is controlled by h-infinity control managed to reduce the effect of the switching transients to the 
higher amount. The red phase of the current waveform varied to only -4A and the blue phase to 4A after the 
switching of the nonlinear load. After that, all three phases settled to their steady-state values with some 
slight variations and very few ripples. When the load was removed from the supply there were very small 
ripples before disappearing and then all the three phases were smooth and stable to the higher amount. 
Moreover, the Loop shaping controller’s response to switching transients was superior as the waveforms of 
the grid current show very slight distortions at the instant of nonlinear load switching, during nonlinear 
loading, and after it was removed. The method has proved to be the fastest in damping switching effects of 
the nonlinear load (see Figure 30). 
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Figure 23. Harmonic spectrum PI Figure 24. Grid current with the PI controller 
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Figure 25. Harmonic spectrum-FLC Figure 26. Grid current with an FLC 
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Figure 27. Harmonic spectrum-H- Figure 28. Grid current with H-infinity control 
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4.3. Discussion 

With the PI controller, the grid current has serious damaging oscillations. The application of the 
FLC has significantly reduced the THD. Additionally, LCL filter, notch filter, and h-infinity control have 
proved to be a better control strategy over PI control and the FLC. Results obtained from simulations show 
the robustness of the robust loop shaping control strategy compared to all the other controllers. 

This experiment has applied the L filter and LCL filter and notch filter and the one in [3] employed 
LC and LCL filters for THD. Both studies have ended up with the reduction of the harmonic currents, the 
LCL filter being the best. The application of a notch filter shows positive effects in [20] as in this research 
although the control strategies are different. Both schemes have managed to significantly reduce the THD of 
the grid current below 5% as required by the IEEE standards. The robustness of LCL in damping harmonics 
can also be seen in [25] as in this study. The application of h-infinity robust control to damp harmonics is 
discussed in [26] similar to this study. In [26], a linear matrix inequality LMI technique is part of h-infinity 
optimization. A robust technique known as linear matrix inequality-linear quadratic regulator (LMI-LQR) is 
employed in [27] to damp voltage harmonics while this study applies robust loop shaping to reduce current 
harmonics. This study and the one in [28] both apply LCL filters with different active damping strategies for 
harmonic reduction. This work uses a notch filter while researchers in [28] prefer washout filter and damping 
coefficient. 


5. CONCLUSION 

An h-infinity robust loop shaping control for controlling the VSI of a PV system in the presence of 
the harmonics induced by nonlinear loads and the inverter is presented in this paper. The proposed 
controller's robustness is demonstrated by the results obtained from experiments with PI, FLC, h-infinity 
control, and robust h-infinity loop shaping control. In the paper, the harmonics in the grid current triggered 
by non-linear loads and the VSI are removed by the LCL and notch filters, as well as the control scheme 
based on robust h-infinity loop shaping. The suitability of the proposed approach has been demonstrated by 
simulation exercises for the two cases. The grid current THD for the two experiments with the proposed 
controller is the lowest, and the grid current waveforms in both cases are the smoothest and most stable for 
the majority of the simulation period as opposed to the other methods. These results prove the robust loop 
shaping method's effectiveness in controlling the VSI. 
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